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Collective flows in heavy-ion collisions from AGS ((2 - ll)A GeV) to SPS ((40, 158)A GeV) ener- 
gies are investigated in a nonequilibrium transport model with nuclear mean- field (MF). Sideward 
ipx), directed v\, and elliptic flows V2 are systematically studied with different assumptions on the 
nuclear equation of state (EOS). We find that momentum dependence in the nuclear MF is im- 
portant for the understanding of the proton collective flows at AGS and SPS energies. Calculated 
results with momentum dependent MF qualitatively reproduce the experimental data of proton 
sideward, directed, and elliptic flows in a incident energy range of (2 — 158)A GeV. 
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I. INTRODUCTION 



Determining the nuclear equation of state (EOS) under 
various conditions has been one of the largest motivations 
of ^I^^Qi^ i"^ )^si(^ in^hes^^dec ades I, .2jJ|j^, .S^^fij^T,, 




the saturation density, EOS gives the bulk properties of 
nuclei such as the binding energy and the radius. While 
the first principle simulations of lattice QCD are possible 
for hot baryon-free nuclear matter and matter at small 
baryon density can be studied by expanding in the power 
series of baryon chemical potential |35j , properties of 
highly compressed matter are still under debate. Thus 
phenomenological studies are necessary to connect the 
experimental heavy-ion collision data with the EOS es- 
pecially for nuclear matter at high baryon densities. In 
high-energy heavy-ion collisions, where nuclear matter in 
a wide range of temperatures and densities are probed, 
many ideas on the EOS and phases have been examined. 
For example, very dense matter is created in recent RHIC 
experiments [s^ suggesting the creation of the gas of de- 
confined quarks and gluons (QGP). In 1970's and 1980's, 
the existence of strong collective flow in heavy-ion col- 
lisions was suggested in hydrodynamics [10, ll| , and it 
was examined in experiments at Bevalac |4|. Since col- 
lective sideward flows are generated in the early stages 
of collisions by the repulsive nucleon potential in nuclear 
matter, the observed strong collective flows were believed 



to signal very large pressure at high baryon densities, 
i.e. hard EOS Q. On the other hand, the real part of 
the nucleon-nucleus potential is already repulsive at the 
normal density at high incident energies, and the role 
of the momentum dependence of nuclear potential on 
the collective flows was extensively studied from around 

1990 USSilk] 

In order to distinguish the momentum 
and density dependence, we need to invoke heavy-ion col- 
lision data in a wide incident energy range. We have now 
systematic collective flow data at various incident ener- 
gies; LBNL BevalaclllElIll, GSI Schwerionen Syn- 
chrotron (SIS) lIlIllllfrMSU NSCLJi|L BNL Alter- 
nating Gradient Synchrotron (AGS) [ll fl l M E 113. l2l| ■ 
CERN Super Proton Synchrotron (SPS) [Hllll^l, and 
BNL RHIC 37]. 

Collective flow data obtained at AGS energies ((2 — 
ll)yl GeV) show a good landmark to determine EOS. 
As demonstrated in Ref. j23, the saturating momen- 
tum dependence of the mean- field (MF), a large num- 
ber of hadronic resonances, and string degrees of free- 
dom are essential in order to explain all of the radial, 
sideward, and elliptic flows at AGS energies. The mo- 
mentum dependence of nuclear potentials in the con- 
text of collective flow was also discussed in Ref. |^ with 
Quantum Molecular Dynamics (QMD) model. It is dis- 
cussed that we can separate the momentum dependence 
by analyzing the so-called balance energy at which the 
flow disappears, and it was later on confirmed in ex- 
periment at NSCL 0] for E;^^ = (55 - 155)A MeV. 
Recently, Danielewicz et al. have also discussed the 
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EOS with these data within Boltzmann equation simula- 
tion 13 HE I27L l28l | , showing that rehable stiffness value 
(^K = 167 — 380 MeV) cannot be uniquely determined 
from currently available collective flow data (f or ^2) up 
to AGS energies {Eir,c = (0.15 - 11)A GeV) On the 
other hand, a comparable description to theirs was ob- 
tained also in Relativistic Boltzmann-Uehling-Uhlenbeck 
(rbuu) by using a Relativistic Mean Field. In rbuu, 
MF is fitted to reproduce the real part of the global opti- 
cal potential in Dirac phenomenology [s^l . In that work, 
a common MF giving K ^ 300 MeV is applied in the en- 
ergy range of (0.25— 11) A GeV. Thus these two works do 
not necessarily provide the same conclusion for the stiff- 
ness. In addition, we still have large ambiguities in the 
MF for hadrons other than nucleons. In order to reduce 
these ambiguities and to pin down the EOS more pre- 
cisely, recently measured flow data at lower SPS energies 
((20 — 80)^ GeV) may be helpful, because higher baryon 
density would be reached at these incident energies. 
Several hadronic transport models, such as RQMD 

M, ELE El, BEM pa mill, RBUU m Isl [32 

ARC j44|. ART |45l|. HSD ^g, UrQMD 0,^3 and JAM 1^ 

have been successfully applied to describe many aspects 
of high energy heavy-ion collisions in a wide range of in- 
cident energies. Transport models without MF effects 
(arc, hsd, jam) can describe bulk observables such as 
transverse mass spectra or rapidity distributions, but 
they cannot explain anisotropic collective flows, which 
are sensitive to MF potentials. Transport models with 
MF effects (rqmd, bem, rbuu, art, urQMo) have been 
successful in explaining anisotropic collective flows in ad- 
dition to bulk observables up to AGS energies. For SPS 
energies, however, the MF effects on collective flows have 
not been seriously investigated yet. 

In this work, we investigate collective flows from 
2 A GeV to 158^ GeV by using a hadronic cas- 
cade model. Jet AA Microscopic Transportation Model 
(jam) combined with a covariant prescription of MF 
(RQMD/S) m. 

This paper is organized as follows. In Sec. m we ex- 
plain our transport model and parameterization of our 
EOS used their. In Sec. CHI we present our results on 
flows on rapidity and transverse distributions as well as 
their excitation functions. In Sec. IIVI we discuss some 
uncertainties in our model. In Sec. [3 we summarize our 
work. 



II. NONEQUILIBRIUM TRANSPORT MODEL 
AND THE EQUATIONS OF STATE 



Heavy-ion collision is a dynamical process of a system 
in which the temperature and density are not uniform 
and the equilibrium is not necessarily reached. Therefore, 
we need dynamical models to describe collisions in order 
to extract static properties of nuclear matter under equi- 
librium. Hydrodynamical description is the most direct 



way to connect the EOS and dynamics. Actually, ideal 
hydrodynamics has succeeded to describe elliptic flow at 
low-p-Ti up to semicentral and around mid-rapidity at 
RHIC '50], (However see Ref. for recent reinterpre- 
tation of the RHIC data.) where the number of produced 
particles is so large that local equilibrium may be easily 
achieved. However the condition of local equilibrium may 
not be satisfied up to SPS energies, and non-equilibrium 
dynamics is required to study the EOS of dense nuclear 
matter through heavy-ion collisions. 

Hadron-string cascade processes are the main source 
of thermalization and particle production up to SPS en- 
ergies. In the increase of incident energy from AGS 
((2-11)^ GeV) to SPS ((20-158)AGeV),the main par- 
ticle production mechanism in hadron-hadron collisions 
evolves from resonance productions to string formations. 
At higher energies, hard partonic interaction (jet pro- 
duction) becomes more important, and the jet produc- 
tion cross section reaches around 20 % of the total cross 
section of pp at RHIC ^3 ■ 

JAM includes all of the above particle and jet pro- 
duction mechanisms, and the applicable incident energy 
range is expected to be enough (for the study of col- 
lective flow, jet production does not matter). Inelastic 
hadron-hadron collisions produce resonances at low ener- 
gies. We explicitly include all established hadronic states 
with masses up to around 2 GeV with explicit isospin 
states as well as their antiparticles, which are made to 
propagate in space-time. At higher energies {^/s > 4 
GeV in BB collisions, y/s > 3 GeV in AIB collisions, 
and -y/s > 2 GeV in AIM collisions), color strings are 
formed and they decay into hadrons after their forma- 
tion time (r ~ 1 fm/c) according to the Lund string 
model PYTHIA 53] . Leading hadrons having constituent 
quarks can scatter within their formation time with other 
hadrons assuming the additive quark cross section which 
is known to be important at SPS energies ji^ . 

It is necessary to include MF effects to explain collec- 
tive flow data, and the MF should have the momentum 
dependence as well as the density dependence in order 
to describe flows in a wide incident energy range. We 
adopt here a simple Skyrme type density dependent MF 
in the zero-range approximation, and a Lorentzian type 
momentum dependent MF j^] which simulates the ex- 
change term (Fock term) of the Yukawa potential. Single 
particle potential U then has the form 
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This MF potential leads to the following total potential 
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TABLE I: Parameter set of density-dependent and momentum-dependent/independent potential. Momentum-dependent hard 
(MH) and soft (MS) potential are taken from Ref. [57l | with simplification (See text for detail). Momentum-independent hard 
(H) and soft (S) potential are taken form Ref. [55|. 
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FIG. 1: Left: Momentum dependence of the single particle potentials Eq. @ for momentum-dependent hard (MH), soft (MS) as 
well as momentum-independent hard (H) and soft (S) are compared with the real part of the global Dirac optical potential ISal . 
Right: Density dependence of total energy per nucleon in Eq. Q for momentum-dependent (MH, MS) and independent (H,S) 
potential. 



energy, through a relation oi U — 6V/df 
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where f(r,p) is the phase space distribution function 



whose integral over p is normalized to the density p{r). 
At zero temperature, the phase space distribution func- 
tion is given as 



/(r,p)= ( -TTp%] p{r) e{pF - \p\) . (3) 



Then the total energy per nucleon is 
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where Fermi momentum is taken to hepp{p) = fi(37r^p/2)^/'^. See ljA22|) for the definition of p in the actual simulations. 
Integrals in Eq. can be obtained analytically as 
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Parameters a, (3 and 7 in Eq. are determined to reproduce the saturation of the total energy per nu- 
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cleon at the normal nuclear density, i.e. E/A\p^pg = 
-16 MeV, and P = p^d[E / A) / dp\p=p^ = MeV/fm^ 
[H^ . The incompressibility is obtained from K = 
'dp'^d'^{E / A)/ dp'^\p^p„. Parameters for hard (H) and soft 
(S) EOS are listed in Table and the density depen- 
dences of the total energy per nucleon are shown in the 
right panel of Fig. ^ 

Ik) 

Parameters Ccx and /i^ are taken to reproduce the real 
part of the global Dirac optical potential (Schrodinger 
equivalent potential) of Hama et al. |38j . in which an- 
gular distribution and polarization quantities in proton- 
nucleus elastic scatterings arc analyzed in the range of 
10 MeV~l GeV in Dirac phenomenology. Single particle 
potential at p = po 
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III. COLLECTIVE FLOWS FROM AGS TO SPS 
ENERGIES 

When two heavy-nuclei collide at high energies at fi- 
nite impact parameters, pressure gradient is anisotropic 
in the initial stages of a collision. As a result, it generates 
anisotropic collective flows. Up to now, several kinds of 
collective flows are proposed to probe high dense mat- 
ter. The first one is the sideward flow (also called di- 
rected flow) {px)i which is defined as the mean value of 
Px^ where x is defined as the impact parameter direction 
on the reaction plane. Sideward flow is mainly gener- 
ated by the participant-spectator interaction. Nucleons 
in the projectile feel repulsive interaction from the target 
nucleus during the contact time of projectile and target. 
This repulsion pushes projectile nucleons out in the posi- 
tive sideward direction if the contact time is long enough. 
When the incident energy is very high, contact time in 
collisions becomes shorter due to the Lorentz contraction, 
therefore sideward flow decreases. At SPS energies, other 
types of collective flows, called as directed (ui) and ellip- 
tic {v2) flows, are mainly measured. These are defined as 
the n-th Fourier coefficient, 



PTdprdydct) 2npTdpTdy 
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is compared to the Schrodinger equivalent potential in 
Ref. [33 in the left panel of Fig. ^ Parameters for 
the momentum dependent potentials are shown as MH 
and MS in Tabled These parameter sets are based on 
Ref. "s?! with simplification in which Coulomb, surface 
and Pauli potentials as well as the zero-point kinetic en- 
ergy of the Gaussian wave packets are dropped because 
their study was focused on nuclear matter below the sat- 
uration density. We have fixed the high energy limit of 
the optical potential, U ^ 77 MeV at i?inc ^ 00, leading 
to a constraint a+f3 = 77 MeV. This constraint generally 
makes EOS stiffer compared to those in Ref. 31] . 

We include the above MF effects into JAM by 
means of simphfied RQMD (RQMD/S) [H frame- 
work. The Relativistic Quantum Molecular Dynamics 
(rqmd) [3^ 0, is a constraint Hamiltonian dy- 
namics, in which potentials are treated in a covariant 
way. RQMD/S 43j uses much simpler and more prac- 
tical time fixation constraints compared to the original 
one [11143,143. For detail, see Appendix A. 

In this work, we take into account potential interac- 
tions only between baryons. Simulation time step size is 
taken to be dt — OA fm/c at all incident energies. We will 
discuss the influence of MF for non-nucleonic baryons on 
the flow analysis and the validity of this treatment in 
section ITVI The violated magnitude of the energy con- 
servation is about 0.4 % in average of time and events. 



where the azimuthal angle (j) is measured from the re- 
action plane. The directed flow vi is the first Fourier 
coefficient of the azimuthal distribution 
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and the elliptic fiow V2 are the second Fourier coefficient 
of the azimuthal distribution 



V2 = (cos2(/)) — 
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These collective flows are reviewed in Ref. '58|. 

The effects of MF in high-energy heavy-ion collisions 
are visible but not very large in single particle spectra, 
such as rapidity distribution dN/dy or transverse mass 
distribution d^N/mTdmj.dy. In this section, we demon- 
strate that MF effects are essential to study anisotropic 
collective flows in the hadron-string transport model JAM 
with MF potentials. 



A. Collective Flows at AGS Energies 

We show proton sideward flow (px) in mid-central 
Au-f Au colhsions at AGS energies C^nc ~ {2 — 11)A 
GeV) together with AGS-E895 data "w] in Fig. Hand in 
the left panel of Fig. 13 We choose the impact parameter 
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FIG. 2: Sideward flows (pa;) of protons in mid-central Au+Au collisions at (2 — 8)A GeV are compared to the AGS-E895 
data [2o|. Lines show the calculated results of Cascade with momentum dependent hard/soft mean-field (MH/MS, left panels), 
Cascade with momentum independent mean-field (H/S, right panels) and Cascade without mean-field (CS). The experimental 
data are shown in both of the left and right panels. 
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FIG. 3: Comparison of calculated sideward fiow (pa;) of protons (left) and pions (right) in llyl GeV Au+Au mid-central 
collisions to AGS-E877 data 1^ . The meaning of the lines is the same as Fig. |2| 



range 4 < 6 < 8 fm in the calculations which roughly cor- 
responds to mid-central collisions in experimental data. 

It is seen that both Cascade and momentum inde- 
pendent soft (S) MF results are inconsistent with the 
data. The magnitude of (pa;) in forward rapidity region 
{y/Upmi — il) is small compared to the data, and the 
slope parameters at mid-rapidity are also smaller than 
that of the data with soft MF. The momentum indepen- 
dent soft MF reduces {px) in forward rapidity region, and 
enhances the slope parameters at mid-rapidity. The for- 
mer is an unfavorable effect in explaining the data, and 
the latter is not enough. With momentum independent 
hard (H) MF, slope parameter is well reproduced, but 
the {px) at forward rapidities are smaller than the data 



especially at i?inc = 2 and ll^GcV. 

Proton sideward flow data are qualitatively reproduced 
with the momentum dependent MF. The momentum de- 
pendent MF pushes up the flow almost linearly as a func- 
tion of rapidity, and it becomes closer to the data, while 
the {px) values at forward rapidities may be a little too 
large compare to the data at E'inc = 4 — 8A GeV. As the 
incident energy increases, MF effects on the slope param- 
eter at mid-rapidity become small, but we can still see 
clear differences at forward rapidities between the results 
with and without momentum dependence. 

Our results suggest the necessity of the momentum de- 
pendence in the MF to give large magnitude emission to 
X direction at forward rapidity. We note that our re- 
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FIG. 4: Transverse momentum de pen dence of the elliptic flow V2 for protons in Au+Au mid-central collisions at (2, 4, 6)A GeV 
are compared to AGS-E895 data [2y| . The meaning of the lines is the same as Fig. |2| 



suits with momentum dependent MF are consistent with 
the previous calculations with MF on the collective flow 
data at AGS energies Is^] as well as SIS ener- 

gies [Ula^- 

The importance of the momentum dependence in the 
MF is also seen in the transverse momentum dependence 
of the proton V2 as shown in Fig. 0] Only if momentum 
dependence is included, we reproduce the strong squeez- 
ing at Sine = GeV. of the pj, dependence. 

In the right panel of Fig. |21 we plot the results of 
sideward flow (p^) for pions in Au-|-Au collisions at 
^-inc = 11^ GeV. The sideward flow (p^) of pions are 
suppressed significantly by momentum dependent MF. 
This may be because pions are trailed by nucleons which 
is affected by MF, giving visible differences. 



B. Directed Flow at SPS energies 

Directed flow vi has been measured at SPS energies 
{Einc = (40, 158)^ GeV) instead of {px) as a function 
of rapidity. In Fig. [S] we compare the rapidity depen- 
dence of proton vi with the data in mid-central Pb-|-Pb 
cohisions at £iji£,= 40^ and 158 A GeV by CERN-NA49 
collaboration |2J| , both of which are deduced by the re- 
action plane method (standard method). One can see 
that momentum dependent MF generally improves the 
description of vi. 

It is interesting to note that the Cascade model over- 
estimates vi for protons contrary to the underestimate of 
(px) at AGS energies. We also see that vi is reduced at 
SPS energies with momentum dependent MF, while (px) 
is enhanced at AGS energies. This is a reverse behavior 
compared to that at lower incident energies. Note also 
that the results with momentum independent MF predict 
larger vi than that of the Cascade results. 

In Fig. |5l the results from the momentum dependent 
MF show a flat behavior at mid-rapidity at 158A GeV. 
The 'wiggle' (a n egative slope of the proton vi near 
mid-rapidity) |29| has been reported at peripheral col- 



lisions 24] . It is interesting to study this in detail in the 
future. 

In Fig. |B| we compare transverse momentum depen- 
dence of vi for protons in Pb-I-Pb collisions at i^inc = 40^ 
and 158^ GeV with the data. We choose rapidity cut 
\y\ < 1.8 for 40A GeV and |y| < 2.1 for 158A GeV ac- 
cording to the experimental cuts. The Pt dependence 
at 158 A GeV is very different from that at AO A GeV. 
Dense baryonic matter is tentatively formed in the cal- 
culations up to around AOA GeV, while many strings are 
formed and hadrons are formed later at 158^ GeV at 
mid-rapidity. As a result, vi does not necessarily grow 
as a function of pr at 158A GeV, because strings do not 
feel MF in our model, and hadrons with large Pt from 
string decay have long formation time in the total CM 
system, and they would have smaller chances to interact 
with other hadrons before strings decay. 

Let us now turn to the pion vi. we show vi for pi- 
ons as a function of rapidity in Fig. |7| and transverse 
momentum in Fig. |S1 at AOA and 158A GeV. It is seen 
that MF effects for pion vi are very small especially at 
mid-rapidities. MF effects are only seen in the forward 
rapidity region for momentum dependent MF. At 40^ 
GeV, in the forward rapidity region, we find reduction 
(enhancement) of vi in the momentum dependent (inde- 
pendent) MF results compared to Cascade ones. This 
comes from the counteraction from protons; momentum 
dependent (independent) MF reduces (enhances) proton 
vi in the mid-rapidity region, and pion vi anti-correlates 
with proton vi. Probably, we need to include pion MF 
for a better understanding of the collective flows at SPS 
energies. 



C. Elliptic Flow at SPS Energies 

Since vi signal becomes small due to the short 
participant-spectator interaction time at high energies, 
the next Fourier coefficient, called as the elliptic flow V2, 
has been discussed more extensively at SPS and RHIC. 
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FIG. 5: Proton directed flows «i as a function of rapidity in mid-central Pb+Pb collisions at Sine ~ 40A GeV (left panel) 
and 158^ GeV (right panel) in comparison with SPS-NA49 data |2^ . Lines show the calculated results of Cascade with 
momentum dependent hard/soft mean-field (MH/MS), Cascade with momentum independent mean-field (H/S) and Cascade 
without mean-field (CS). 
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FIG. 6: Proton directed flows vi as a function of transverse momentum in mid-central Pb+Pb collisions at _Einc = 40A GeV 
(left panel) and 158 A GeV (right panel) are compared with SPS-NA49 data [23 |. The meaning of the lines is the same as Fig.|3 



At these energies, the participants form an almond-like 
shape in the transverse plane after the spectators go 
through, and this almond shaped participants start to 
expand more strongly in the x (shorter axis of the al- 
mond) direction, due to the higher pressure gradient if 
the participants are well thermalized. This expansion is 
known to lead to the enhancement of in-plane particle 
emission, i.e. positive elliptic flow V2- 

In Figs. El and 1101 we plot the results of the rapid- 
ity and transverse momentum dependence of V2 for pro- 
tons at SPS energies (40A and 158A GeV) together with 
the SPS-NA49 data At SPS energies, the Cascade 
model generally explains the proton V2 data qualitatively, 
including the flat behavior of ^2(2/) at mid-rapidities at 



158 A GeV and the approximate linear pj- dependence of 
V2iPT)- One exception is the missing collapse of V2{y) at 
mid-rapidity at 40 A GeV. This collapse seen in the NA49 
data may be an indication of a first order phase transition 
at high baryon densities achieved in the Pb-|-Pb collisions 
at iOA GeV (H- 

Effects of MF are small for proton V2 at SPS energies. 
Elliptic flow is most easily generated in the early stages 
of the collisions, since spatial anisotropy is the largest. 
However, at SPS energies, string excitations dominate 
particle production at early times in the model and those 
strings are not affected by the nuclear mean field. That 
is the reason why MF effects are small at SPS energies 
in our results. 
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FIG. 7: Pion directed flows vi as a function of rapidity in mid-central Pb+Pb collisions at Einc = 40A GeV (left panel) and 
158A GeV (right panel) are compared with SPS-NA49 data |2^. The meaning of the lines is the same as Fig.|5] 
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FIG. 8; Pion directed flows vi as a function of transverse momentum for < y < 1.5 in Pb+Pb collisions at _Einc ~ ^OA 
GeV (left panel) and 158^4 GeV (right panel) are compared with SPS-NA49 data i24il . The meaning of the lines is the same as 
Fig.ini 



Rapidity and transverse momentum dependence of the 
pion V2 are shown in Figs.lllland ll2l respectively. Rapid- 
ity dependence at 158 A GeV and transverse momentum 
dependence at low Pt {pt < 1 GeV/c) at 40A and 158A 
GeV are well explained by the Cascade model as well as 
by the momentum dependent /independent MF models, 
and we do not find any significant MF effects for these ob- 
servables. By contrast, we do not see the collapse of V2{y) 
at mid-rapidities seen in the iOA GeV NA49 data, and 
we underestimate V2 at high p^.. The former corresponds 
to the collapse of proton V2{y) mentioned before. Mo- 
mentum independent MF enhances pion V2{y) slightly, 
but this is in the reverse direction to explain the data 
at 40 A GeV. The strong increase of V2{pt) up to around 
Pt ^ 2 GeV/c is also seen at RHIC energies, and this 



behavior is discussed as an indication of hydrodynamical 
evolution [5(tI |. 

We now turn to the discussion of the difference be- 
tween vi and V2- We have shown that MF effects on vi 
are rather strong, but V2 is relatively insensitive to MF 
at SPS energies. This may come from the difference of 
developing time between vi and V2 ■ The directed flow vi 
at mid-rapidities is mainly generated by the interaction 
between participants and spectators in the early stage of 
the collision, where baryon density is the highest. On 
the other hand, V2 in our model is generated in the late 
stage until the time reaches the order of nuclear radius, 
where densities are not very high, but it is not formed in 
the early stage in our model. This is because our current 
hadronic transport approach does not have large partic- 
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FIG. 9: Proton elliptic flows V2 as a function of rapidity in mid-central Pb+Pb collisions at 40^4 GeV (left) and 158^4 GeV 
(right). Lines show the calculated results of Cascade with momentum dependent hard/soft mean-field (MH/MS), Cascade 
with momentum independent mean-field (H/S) and Cascade without mean-field (CS). Experimental data are taken from SPS- 
NA49 pl|. 




FIG. 10: Proton Elliptic fiows V2 as a function of transverse momentum. The meaning of the lines is the same as Fig.|5| 



ipant pressure in the early stages of the collisions, as we 
do not explicitly include MF for strings and partonic in- 
teractions. In a hydrodynamic picture, V2 develops from 
very early times due to thermal pressure. This is a strik- 
ing difference between our approach and hydrodynamics 
as previously studied in Ref. HJ] . 



D. Elliptic flow excitation functions from AGS to 
SPS Energies 

When the incident energy is not high enough, specta- 
tors squeeze participants out of the reaction plane due 
to the repulsive nuclear interactions at 0.2 < Ei^c ^ 4:A 
GeV. This squeezing leads to a negative value of the ellip- 
tic flow of nucleons {v2 < 0). The elliptic flow, therefore, 
shows the strength of the repulsive interaction at lower 



energies. On the other hand, elliptic flow becomes posi- 
tive at higher energies, because there is no such squeezing 
effect due to the Lorentz contraction. Elliptic flow gives 
a information how much pressure is generated at higher 
energies. 

In Fig. ll3[ we show the incident energy dependence of 
proton V2 in mid-central collisions with measured data 
(-0.1 < y < 0.1 for AGS, < y < 2.1(0 < y < 
1.8) for SPS 158A (40A) GeV) [H H^. Rapidity cut 
\y\ < 0.2ypi-oj has been used in calculations. Experimen- 
tal data clearly show the evolution from squeezing to al- 
mond shaped participant dynamics. With both Cascade 
and momentum independent soft MF (S), we cannot ex- 
plain strong squeezing effects at lower energies. The cal- 
culated V2 values for momentum independent MF (H,S) 
and Cascade are generally larger than data at AGS en- 
ergies. Momentum dependent MF (MII,MS), which is 
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FIG. 11: Pion elliptic flows V2 as a function of rapidity. The meaning of the lines is the same as Fig. El 
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FIG. 12: Pion elliptic flow V2 as a function of transverse momentum. The meaning of the lines is the same as Fig.|^ 



repulsive in the incident energy range under considera- 
tion, pushes down the eUiptic flow significantly. We qual- 
itatively reproduce the incident energy dependence from 
AGS 19] to SPS 13 energies. 

Calculated results with both MH and MS are smooth 
as a function of beam energy, while the data at Ei^c = 
AO A GeV has a dip Confirmation of data is nec- 

essary to examine the incident energy dependence of V2, 
whether it is a monotonic function or has a dip at around 
Einc ^ 40^ GeV by looking at the missing data points. 

In our results with momentum dependent potentials, 
the stiffness dependence of V2 is smaller than that in the 
Boltzmann Equation Model (bem) [Ulllllil. In the 
RQMD/S framework with the relativistic distance ffj, 
the interaction between the projectile and target nucle- 
ons are suppressed at high energies by the factors mi/pi 
and rrij/p^ in the potential derivatives in Eqs. (|A25I) and 
(jA26|) . For momentum dependent potentials, we have 
the relative momentum vector pij in Eq. (|A29|) . which 
can compensate the suppression factor in Eq. (jA26|) . For 



momentum independent potentials, on the other hand, 
the pair velocity Pij in Eq. (|A20|) is very small for nu- 
cleon pairs between the projectile and target, and there 
is no enhancement factor to compensate the above sup- 
pression in the derivatives of the relativistic distance ffj 
in Eqs. l|A27p and ljA28|l . This suppression does not hap- 
pen in BEM, and t hey find significant stiffness dependence 
in Refs. 0, |23, while we do not see strong stiff- 
ness dependence. In the case of momentum independent 
potentials, our results are closer to the Cascade results 
compared to those in Ref. ji^, I3| ■ This difference also 
comes from the above suppression between the projectile 
and target nucleons. Essential reason for those difference 
is that in rqmd or RQMD/S, potentials are regarded as 
Lorentz scalar. Possible other model dependence will be 
discussed in the next section. 
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FIG. 13: Incident energy dependence of proton elliptic flow at mid-rapidities in mid-central heavy-ion collisions from lA GeV 
to 158A GeV. Dotted, bold-solid, dashed, and dot-dashed and thin-solid lines show the results of Cascade with MH, MS, H, S 
mean-field, and Cascade without mean-field (CS), respectively. The experimental data are taken from LBL-EOS, AGS-E895, 
E877 from Ref. (igj and SPS-NA49 from Ref. 



IV. MODEL UNCERTAINTIES 

In the previous section, it has been shown that the 
momentum dependent hard or soft MF improved the de- 
scription of the collective flow data from AGS to SPS 
energies. However, there are some uncertainties in our 
calculations for the study of collective flows. 

First, let us consider the effects of the MF for non- 
nucleonic baryons. Strange baryons, resonance hadrons 
or anti-baryons are expected to feel MF, which may be 
different from that for nucleons. In the previous section, 
we have assumed that all the baryons feel the same MF, 
and this treatment would give a rough estimate of a max- 
imum MF effects, since, for example, MF for A's or A's 
is generally expected to be smaller than that of nucleons. 
On the other hand, if we include MF only for nucleons, 
we may get a rough estimate of a minium baryonic MF 
effects. 

In Figs. El and ^1 we compare the results with and 
without MF for non-nucleonic baryons. We can read 
in the left panel of Fig. E| that ignoring MF for non- 
nucleonic baryons (specified as "N" in the figure) at 2A 
GeV reduces both the (px) slope and the strength at the 
forward rapidities by about 20 % compared to the case 
of "B" in which all baryons feel MF. But slope remains 
the same at llA GeV among "N" and "B". On the other 
hand, it is seen in the right panel of Fig.^jthat MF only 
for nucleons is not enough to suppress vi at SPS energies. 
One can also see some differences of V2 in "B" and "N" 
in Fig. 1151 up to AGS energies. The experimental data 



lie between "B" and "N" except for 40 A GeV, suggest- 
ing that MF for excited baryons are smaller than that of 
nucleons. 

Next, we have checked the time step size dt depen- 
dence. Since update of MF after each collision requires 
a huge calculation time, we evaluate MF only at each 
time slice. When a baryon collides in one time step, that 
baryon is propagated with MF until the collision time and 
we ignore the MF after the collision before it is formed. 
In the time step of baryon formation, displacements by 
the MF for pi and are evaluated by using the MF at the 
next time slice. This treatment is valid up to the first or- 
der in dt when one baryon collides and/or is formed once 
in one time step. In the later stages, this prescription 
is expected to work well because of the low collision fre- 
quency. In the early stages, many collisions make strings 
and resonances which do not feel MF, so our prescrip- 
tion may not be too bad. In the middle stages, however, 
it may be possible that elastic scatterings are frequent 
enough and baryons keep to feel MF after collisions in 
each time step. Thus we need to analyze the collision 
frequency effects on constructing flows by reducing dt. 
In Figs. 1141 and llSI we plot the results with different time 
step sizes dt— 0.05, 0.1 and 0.5 fm/c. For {p^} at AGS en- 
ergies (left panel of Fig. ll4|) . all the results with different 
time step size agree well with each other, and we can- 
not distinguish these lines. Time step size dependence 
of the V2 as shown in Fig. E] still gives us the confirma- 
tion of the convergence of the numerical results. For vi 
at SPS energies (right panel of Fig. I14|l . only very small 
differences can be seen between the results with dt = 0.5 
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FIG. 14; The sideward (directed) flows of AGS (ieft) and SPS (rigfit) energies are compared witli diflerent assumptions for 
tfie mean-fields. Momentum dependent soft (MS) type mean-field has been used in these calculations. The first term in the 
parenthesis means that all baryon (B) or only nucleon (N) are included for the MF. 

The second term (0.05/0.1/0.5) denotes simulation time step size dt (fm/c). 




FIG. 15: The elliptic flows are compared with different assumptions for the mean-fields. Meaning of the lines is the same as 
Fig.m 



fm/c and dt < 0.1 fm/c. We conclude that dt — 0.1 fm/c 
which has been used as a default throughout this work is 
small enough to perform a reliable calculations. 

Finally, we would like to address the problem of the un- 
certainties of the transport model itself. In addition to 
the ambiguities in introducing collision terms, the equa- 
tions of motion depend on the model treatment. It is 
not trivial at all to construct equations of motion of rela- 
tivistic particles during heavy-ion collisions based on the 
potential or the MF giving an appropriate EOS. At rel- 
ativistic energies, there are several ways proposed so far 
to introduce the potential effects. 



1. Relativistic Mean Field (RMF) having Lorentz 
scalar Us and vector terms (rbuu [H, Is^). 
The scalar and the vector time-component are eval- 
uated in the local rest frame, and by the Lorentz 
transformation we can get C/^ in the calculation 
frame. Thus neglecting a nonlocality in time, this 
evaluation of the potential is practically covariant. 
In this approach, however, we need to introduce 
strong cut-off for the coupling of vector meson and 
baryons [S^, since the vector potential effects lin- 
early increase as a function of incident energy. 
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2. Lorentz scalar re-interpretation of non-relativistic 
potentials (buu [sJI, BEM [23|). In the BUU 
model the Lorentz scalar MF Us is obtained 
from the non-relativistic MF U in the local rest 
frame through the relation, 

e{p,p) = v^(m + Us{p, p)Y 

= ^Jm?+p^ + U{p,p) , (10) 

where p is the baryon density [sj- For the mo- 
mentum independent MF in BEM in Refs. p ^ l28 | . 
the scalar potential is directly given so as to fit the 
EOS, and the scalar density is used for p. They do 
not have any vector terms increasing at high en- 
ergies, and the potential effects become mild com- 
pared to the RMF treatment. For example, the 
derivative of the above single particle energy gives 
rise to the factor {m + Us)/e \n front of the Ug 
derivative, and suppresses the potential effects. 

3. Combination of the Lorentz scalar and non- 
relativistic type density dependent potentials 
(bem 2d|, I23, I23I). This approach is adopted in 
Ref. 1261 and in the momentum dependent MF in 
Refs. [23, ■ The single particle energy is given as 

e{p,p)=m+f dp'v*{p',p) + U{p) , (11) 
Jo 

v*{p,p) = ^= P (12) 
yp2 ^ (m*{p,p) f 

The derivative of U (p) does not come with the sup- 
pression factor such as m/e, and generates strong 
effects at high energies, where both of the density 
and density derivative become large. 

4. Constraint Hamiltonian dynamics (rqmd j3^ . 
RQMD/S "4^). In rqmd and RQMD/S, particle 
velocity and force are not given by the derivatives 
of the single particle energy, but by the derivatives 
of the total Hamiltonian, as shown in the Appendix 
in the case of RQMD/S. Thus the relation to other 
MF models described above is not straightforward. 
However, the potential Vi in the on-mass-shell con- 
straint (|A8I) is introduced as Lorentz scalar, and 
we have suppression factor m/p^ in the equations 
of motion, (|A14p and l|A15|l . These observations 
suggest that the rqmd and RQMD/S would give 
results similar to those in Lorentz scalar MF mod- 
els, such as BUU [33. Another difference from other 
MF models exists in the nuclear density profile. 
One nucleon is represented by one Gaussian packet 
rather than many test particles, then the nuclear 
diffuseness becomes generally larger in QMD-type 
models. This may generate artificial surface ef- 
fects at large impact parameters or in light-ion col- 
lisions. However, central and mid-central collisions 
of heavy-nuclei are expected to be well described. 



as in the case of various cascade models, in which 
one particle is used for one hadron. 

We are not very sure which is the best way to include 
the potential effects in high-energy heavy-ion collisions. 
Further formal developments on transport models in re- 
lation to nuclear EOS would be necessary, and at the 
same time, phenomenological studies of heavy-ion colli- 
sions are required to verify the validity of models and to 
elucidate the EOS. From the latter point of view, system- 
atic study in a wide range of incident energy is needed, 
since the above uncertainties are closely related to the 
Lorentz transformation properties, whose effects would 
vary drastically as the incident energy varies. The inci- 
dent energy range from AGS to SPS energies studied in 
this work may provide a good benchmark test for trans- 
port models and the EOS. 

V. SUMMARY 

We have investigated collective flows in heavy-ion colli- 
sions from AGS ((2-ll)AGeV) to SPS ((40, 158)A GeV) 
energies by using a combined framework of hadron-string 
cascade (jam) 49] and covariant constraint Hamiltonian 
dynamics (RQMD/S) I'^ JAM, various particle pro- 
duction mechanisms are taken into account — production 
and decay of resonances and strings, jet production and 
its fra gme ntation. Momentum dependence of the MF is 
fitted [53 to the real part of the Schrodinger equivalent 
global optical potential of Hama et al. [s^ in a Lorentzian 
form Eq.JQ). Saturation properties are fitted by intro- 
ducing the density dependent potential of Skyrme-type 
in the power series of p. Calculated results of Cascade, 
Cascade with momentum dependent MF, and Cascade 
with momentum independent MF are compared with the 
data of sideward (px), directed wi, and elliptic V2 flows as 
a function of rapidity, transverse momentum and beam 
energy from AGS to SPS. Generally, results with momen- 
tum dependent MF reasonably well explain the trend of 
the data for proton flows. We note that it is for the first 
time to explain anisotropic proton collective flow data 
of heavy-ion collisions from AGS to SPS in one frame- 
work consistently. Without momentum dependence in 
MF, we cannot reproduce the strong enhancement of the 
sideward flow at E^inc = (2— 11) A GeV, strong squeezing 
seen in V2 for E^inc ^ 4A GeV, and the suppression of 
proton vi at Einc = 40A and 158A GeV. 

Our new model — hadron-string cascade with momen- 
tum dependent MF — provides an improved description 
for collective flows in mid-central collisions from AGS to 
SPS energies. The present analysis implies that the ef- 
fects of the momentum dependent potential is large up 
to the SPS energies. 

There are still many problems to pin down the equa- 
tion of state of dense nuclear matter from heavy-ion data. 
First, we have made an assumption that MF is taken into 
account only for baryons and all the baryons feel the same 
MF. It is interesting to extend the present work to discuss 
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the MF effects for mesons and different MF for hyperons 
and resonance hadrons, and fook at the A or kaon flow 
data. Secondly, we cannot make soft EOS {K ^ 200 
MeV) in the present form of MF which is consistent with 
the optical potential. When the momentum dependence 
is fitted to the optical potential by Hama et al. in the 
Lorcntzian form, the EOS necessarily becomes relatively 
stiff in combination with the Skyrme-type density depen- 
dent form as shown in Table |l| The small sensitivity on 
EOS with momentum dependence appeared in this work 
may be suggesting that the probed EOS range is not wide 
enough. Finally, the model dependence of the MF treat- 
ment has to be cleared in order to obtain model indepen- 
dent EOS information. For this purpose, it is necessary 
to test various MF treatment in one framework. In the 
present model, we need to modify the on-mass-shell con- 
straint to include Lorentz vector potentials or potentials 
of other types. It can be a breakthrough for a transport 
model-independent discussion of EOS. 
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APPENDIX A: RQMD/S FORMALISM 

Here we briefly summarize the ROMP / S formalism de- 
veloped by Maruyama et al. in Ref. |43 for completeness. 
The original RQMD formalism is initiated by Sorge et 
al. in Ref. [s^]- RQMD(/S) is based on the constrained 
Hamiltonian dynamics 59] which is formulated in a man- 
ifestly covariant way. We use 4-vectors and for 
the description of the N particle system. Therefore, we 
need to have 2 A'' constrains (/)i(i — 1, . . . ,2A^) as physi- 
cal phase space is 6N dimension. Now our Hamiltonian 
may be constructed from the constraints 0, and Lagrange 
multiplier Ui from the Dirac's constraint Hamiltonian for- 
malism 



2N-1 



(Al) 



1=1 

The equations of motion are then 

2N-1 



dr 



2N-1 „ 



(A2) 
(A3) 



where the Poisson brackets are defined as 

dA dB dA dB^ 



{A,B}^Y.{ 



(A4) 



(A5) 



i,j,k = 1, . . . ,iV , ^, J/ = 0,1,2,3 



and the sign means the weak equality initiated by 
Dirac When we require that constraints (f>i should 
conserve in time, then they fulfill 



dr 



dr 



2N-1 

^ Uj{(j)i,(j)j} « 0. 



(A6) 



Particle trajectories in QN phase space is uniquely de- 
termined by the equations of motion Eqs. (|A2|) and (|A3|) 
together with Eqs. ljA10|) and ljAll|) when 2N constraints 
are given. 

We use the following 2N constraints in RQMD/S 




z = l,. 

i = N ■ 



.N 



1,. 



,2A^. 



(A7) 



First A^ constraints are the on-mass-shell constraints 



2m,Vi w . 



(A8) 



Remaining N conditions constrain the time fixation of 
the particles. In original RQMD time fixation ^39, 40, 4^, 
A^ X A^-matrix has to be solved numerically at each time 
step to deduce inverse matrices. Moreover, if particle 
production or annihilation occurs, the time fixation is vi- 
olated and a initial qi of produced particles satisfying the 
constraints and energy conservation has to be imposed. 

Maruyama et al. introduced simplified time fixation in 
RQMD/S with global time parameter r in Ref. ^3 as 



Xt = d ■ {qi - qM) ^ , 
Xn = d ■ qN - T KiO^ 



(A9) 



where, a is a 4-component vector corresponding to (1, 0) 
at the rest frame of the particle, and qi is space-time 
coordinate of the i-th particle. The constraints Eq. (|A9|) 
is able to be kept in the case of particle productions. 

Since the constraints (j)i{i = 1, . . . , 2A^ — 1) does not 
depend explicitly on r, Lagrange multiplier Ui(r) can be 
solved as 



02 N 



dr 



2N,i, 



i = 1, . . . ,2iV- 1, 



where C^/ = «, j = 1, ■ 



,2N. 



(AlO) 
(All) 



The matrix C (inverse of matrix C^^) must exist, be- 
cause we only allow for the r-dependent 2Nth constraint 
functions which combine with 2A'^ — 1 constraints. Fur- 
thermore, C can be obtained analytically, if we replace 
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in the potential Vi with the kinetic energy -^/p^ + m| . 
This is a great advantage in the point of view of CPU 
time. One obtains the RQMD/S Hamiltonian 



N 



where 



1 



2m,: 



The equations of motion are then 



dr 

dpi 
dr 



dH 
dpi 



Pj_ 



rrij dVj 



dH 
dr. 



f^i P° 9p^ 

ruj dVj 



(A12) 
(A13) 

(A14) 
(A15) 



In actual calculations, we have replaced with 
the kinetic ener gy y/ pf + mf in the denominators of 
Eqs. (|A14|1 and ljA15|) after evaluating all the derivative 
terms for simplicity. This approximation would be valid 
in the relativistic energy region, where the kinetic energy 
is much larger than the potential Vi. 

Relative distance — Vi — Vj and pij — Pi — Pj in 
the potentials should be replaced by the squared 4- vector 
distance with Lorentz scalar as. 



-iTij 



p. .2 



{pi: 



P 2 



(A16) 
(A17) 



where pij = Pi - Pj , qij = qi - qj 
We note that in non-relativistic limit, 



Pij = Pi 

2 

-qTij 



+ Pj. 

■> r . 



This assumption takes into account the contraction of 
longitudinal direction and we can avoid unphysical com- 
pression. In the actual simulations, we use the following 
expression 



2 



-PTij 



— IJ 



^2 



0\2 



(A18) 

n2 \ 2 



Pt+P' 



(A19) 



where, the velocity and 7-f actor between i- and j-th par- 
ticle are given by 



ft. 



Pi + Pj 

P°+p', ' 



1 



1^3 



(A20) 



We now write down the explicit form of the equations 
of motion in RQMD/S which is used in the actual sim- 
ulation. As explained in Sec. ^1 we use the following 



potentials 

V — y^(Vsky i + Vmom i) 

i 



E 



2po 



(I+7K 

2po ^ TTfe 



{Pi}' 



PrjlPkf 



(A21) 



where [pi) is obtained from a convolution of Gaussian 
wave packet: 



{Pi) = E / drpi{r)pj{r) = ^ 



Pi] 



E 



(47ri)^ 



r2. 

■expi-^ 



(A22) 



The width parameters L =2.05(MHV2.1(MS), 1.08(H 
and S) fm^ are taken from Refs. |55Ll57l |. The equations 
of motion Eqs. ljA14p and (|A15|I then become 



dr p\ 

dpi 
dr 



Ea 



dpi 



^p^ 



(A23) 
(A24) 



where. 



Dij ^\-^]P^3 



1 + 7 Po 1 



7-1 







''rrii 


rrij \ 


( 4L, 




k'W' 





(fc) 



fc=l,2 



[Ptj/Pk 



(A25) 



1 f mi 



fc=l,2 



r(k) 



p.k^ J 1+ [Pij/flk]^' 

(A26) 



The result of the differentials are [6l| 



(9f? 
dpi 



nirrj ■ ftj)<j + 7»j-(^'»J • fti) ( ftj - ^ 



(A27) 



dr^ 

9pI 
dpt 



27; 



,0 



(A28) 



(A29) 
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Finally, let us check the non-relativistic limit to con- 
firm the validity of Eqs. 1A12II - HA15I) . We define kinetic 
energy as £i = ~ rriiC^ (here we write the speed of 
light c explicitly), Indeed Hamiltonian Eq. (|A12|) and the 
equations of motion have the correct non-relativistic limit 
as 

^1 1 , \ 

E(£:'-^-^^)-^-^N.R. , (A30) 



c — '■OO 



dri 9i/N.R. Pi , dV. 



dr dpi rrii 9p 
dpi _ dH^.n. ^dVj 



W. Scheid, H. Miiller and W. Greiner, Phys. Rev. Lett. 
32, 741 (1974). 

J. Kapusta and D. Strottman, Phys. Lett. B 106, 33 
(1981). 

H. Stocker et al, Phys. Rev. C 25, 1873 (1982). 

H. A. Gustafsson et al, Phys. Rev. Lett. 52, 1590 (1984). 

G. F. Bertsch and S. Das Gupta, Phys. Rep. 160, 189 

(1988). 

C. Gale, G. Bertsch, S. Das Gupta, Phys. Rev. C 35, 
1666 (1987). 

G. M. Welke, M. Prakash, T. T. S. Kuo, S. Das Gupta, 
and C. Gale, Phys. Rev. C 38, 2101 (1988). 
L. P. Csernai, G. Fai, C. Gale, and E. Osnes, Phys. Rev. 
C 46, 736 (1992). 

S. Soff, S. A. Bass, C. Hartnack, H. Stocker and 

W. Greiner, Phys. Rev. C 51, 3320 (1995). 

J. Chance et al. (EOS Collab.), Phys. Rev. Lett. 78, 2535 

(1997). 

J. C. Kintner et al., Phys. Rev. Lett. 78, 4165 (1997). 
K. G. R. Doss et al, Phys. Rev. Lett. 57, 302 (1986). 
N. Bastid et al. (FOPI Collab.), Nucl. Phys. A622, 573 
(1997). 

A. Andronic et al. (FOPI Collab.), Phys. Rev. C 67, 
034907 (2003). 

G. Stoicea et al. (FOPI Collab.) and P. Danielewicz, 
Phys. Rev. Lett. 92, 072303 (2004). 

R. Pak et al., Phys. Rev. C 54, 2457 (1996). 

P. Chung et al. (E895 Collab.) and P. Danielewicz, Phys. 

Rev. C 66, 021901(R) (2002). 

J. Barrette et al. (E877 Collab.), Phys. Rev. C 55, 1420 
(1997); 56, 3254 (1997). 

C. Pinkenburg et al. (E895 Collab.), Phys. Rev. Lett. 83, 
1295 (1999). 

H. Liu et al. (E895 Collab.), Phys. Rev. Lett. 84, 5488 
(2000). 

J. L. Klay et al. (E895 Collab.), Phys. Rev. Lett. 88, 
102301 (2002). 

H. Appelshauser et al. (NA49 Collab.), Phys. Rev. Lett. 
80, 4136 (1998). 



[23] H. Appelshauser et al. (NA49 Collab.), Phys. Rev. Lett. 

82, 2471 (1999). 
[24] C. Ah et al. (NA49 Collab.), Phys. Rev. C 68, 034903 

(2003). 

[25] P. K. Sahu, W. Gassing, U. Mosel and A. Ohnishi, Nucl. 

Phys. A672, 376 (2000). 
[26] P. Danielewicz, et al., Phys. Rev. Lett. 81, 2438 (1998). 
[27] P. Danielewicz, Nucl. Phys. A673, 375 (2000). 
[28] P. Danielewicz, R. Lacey, W. G. Lynch, Science 298, 

1592 (2002). 

[29] R. J. M. Snellings, H. Sorge, S. A. Voloshin, F. Q. Wang, 
and N. Xu, Phys. Rev. Lett. 84, 2803 (2000). 

[30] D. Persram and C. Gale, Phys. Rev. C 65, 064611 (2002). 

[31] A. B. Larionov, W. Gassing, C. Greiner and U. Mosel, 
Phys. Rev. C 62, 064611 (2000). 

[32] T. Maruyama, W. Gassing, U. Mosel, S. Teis and K. We- 
ber, Nucl. Phys. A573, 653 (1994). 

[33] P. K. Sahu and W. Gassing, Nucl. Phys. A712, 357 
(2002). 

[34] H. Stocker, E. L. Bratkovskaya, M. Bleicher, S. Soff and 
X. Zhu, J. Phys. G 31, S929 (2005). 

[35] See, for example, K. Redlich, F. Karsch and A. Taw- 
fik, J. Phys. G 30, S1271 (2004) arXiv:nucl-th/0404009 ; 
F. Karsch, Prog. Theor. Phys. Suppl. 153, 106 (2004) 
arXiv:hep-lat/0401031 ; C.R. AUton , M. Doring, S. 
Ejiri, S.J. Hands, O. Kaczmarek, F. Karsch, E. Laer- 
mann, K. Redlich, Phys. Rev. D 71, 054508 (2005) 
arXiv:hep-lat/0501030 , and references therein. 

[36] See, for example, Nucl. Phys. A698 (2002), and refer- 
ences therein. 

[37] C. Adler et al. (STAR Collab.), Phys. Rev. C 66, 034904 
(2002); B. B. Back et al. (PHOBOS Collab.), Phys. Rev. 
Lett. 89, 222301 (2002); S. S. Adler et al. (PHENIX Col- 
lab.), Phys. Rev. Lett. 91, 182301 (2003); J. Adams et 
al. (STAR Collab.), Phys. Rev. Lett. 92, 062301 (2004). 

[38] S. Hama, B. C. Clark, E. D. Cooper, H. S. Sherif, and 
R. L. Mercer, Phys. Rev. C 41, 2737 (1990). 

[39] H. Sorge, H. Stocker and W. Greiner, Ann. of Phys. 192, 
266 (1989). 



17 



[40] H. Sorge, Phys. Rev. C 52, 3291 (1995). 

[41] H. Sorge, Phys. Lett. B 402, 251 (1997); Phys. Rev. Lett. 

82, 2048 (1999). 
[42] T. Maruyama et al, NucL Phys. A534, 720 (1991). 
[43] T. Maruyama et al, Prog. Theor. Phys. 96, 263 (1996). 
[44] Y. Pang, T. J. Schlagel and S. H. Kahana Nucl. Phys. 

A544, 435c (1992). 
[45] B. A. Li and C. M. Ko, Phys. Rev. C 52, 2037 (1995). 
[46] W. Ehehalt and W. Gassing, Nucl. Phys. A602, 449 

(1996). 

[47] L. A. Winckelmann et al, Nucl. Phys. A610, 116c 
(1996). 

[48] S. A. Bass et al, Prog. Part. Nucl. Phys. 41, 255 (1998). 

[49] Y. Nara, N. Otuka, A. Ohnishi, K. Niita, and S. Chiba, 
Phys. Rev. C 61, 024901 (2000). 

[50] P. F. Kolb, P. Huovinen, U. W. Heinz and H. Heiselberg, 
Phys. Lett. B 500, 232 (2001); P. Huovinen, P. F. Kolb, 
U. W. Heinz, P. V. Ruuskanen and S. A. Voloshin, 
Phys. Lett. B 503, 58 (2001); P. F. Kolb, U. W. Heinz, 
P. Huovinen, K. J. Eskola and K. Tuominen, Nucl. Phys. 



A 696, 197 (2001); P. Huovinen, arXiv:nucl-th/05 05036| 
T. Hirano, Phys. Rev. C 65, 011901 (2002). 
[51] T. Hirano and M. Gyulassy, arXiv:nucl-th/0506049 
[52] X. -N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 
(1991). 

[53] T. Sjostrand et al, Gomp. Phys. Gomm. 135, 238 (2001). 
[54] S. Soff, S. A. Bass, M. Bleicher, H. Stcker, and 

W. Greiner, arXiv:nucl-th/9903061 
[55] J. Aichelin and H. Stocker, Phys. Lett. B 176, 14 (1986). 
[56] J. Aichelin, A. Rosenhauer, G. Peilert, H. Stoecker, and 

W. Greiner, Phys. Rev. Lett. 58, 1926 (1987). 
[57] T. Maruyama et al, Phys. Rev. G 57, 655 (1998). 
[58] J. -Y. OUitrault, Nucl. Phys. A638, 195c (1998). 
[59] A. Komar, Phys. Rev. D 18, 1881; ibtd. 1887; tbid. 3617 

(1978). 

[60] P. A. M. Dirac, Rev. of Mod. Phys. 21, 392 (1949). 
[61] Y. Hirata et al, Nucl. Phys. A707, 193 (2002); Y. Hirata, 
Ph.D. thesis, Hokkaido University (2000). 



